arXiv:1508.02264v2 [quant-ph] 4 Jan 2016 


Generation of cluster states in optomechanical quantum systems 

Oussama Houhou, 1, 2 , E] Habib Aissaoui, 1 ’!]] and Alessandro Ferraro 2 ’@ 
1 Laboratoire de Physique Mathematique et Subatomique (LPMS), 

Universite de Constantine 1, Constantine, Algeria 
2 Centre for Theoretical Atomic, Molecular and Optical Physics, 

School of Mathematics and Physics, Queen’s University, Belfast BT7 INN, United Kingdom 

(Dated: August 05, 2015) 

We consider an optomechanical quantum system composed of a single cavity mode interacting 
with N mechanical resonators. We propose a scheme for generating continuous-variable graph states 
of arbitrary size and shape, including the so-called cluster states for universal quantum computation. 
The main feature of this scheme is that, differently from previous approaches, the graph states are 
hosted in the mechanical degrees of freedom rather than in the radiative ones. Specifically, via a 2 N- 
tone drive, we engineer a linear Hamiltonian which is instrumental to dissipatively drive the system 
to the desired target state. The robustness of this scheme is assessed against finite interaction times 
and mechanical noise, confirming it as a valuable approach towards quantum state engineering for 
continuous-variable computation in a solid-state platform. 
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I. INTRODUCTION 

Recent experimental advances have shown that vari¬ 
ous types of mechanical oscillators can operate deeply 
in the quantum regime EH, promoting these systems 
to interesting candidates for quantum technologies. In 
particular, different cooling techniques have succeeded in 
bringing these oscillators close to their ground state [U- 
|9j, whereas the ability to realize a coherent radiation- 
pressure interaction between electromagnetic and me¬ 
chanical degrees of freedom has allowed for the realization 
of genuine quantum processes, such as quantum state 
transfer Ea and the generation of squeezing EMI and 
entanglement El- These achievements, together with 
the possibility to scale up the number of involved oscil¬ 
lators 154181] . pave the way for more advanced quan¬ 
tum information applications such as the engineering of 
quantum dissipation El, quantum many-body sirnula- 
4221 , and quantum information processing in gen- 


In this context, Schimidt et al. |24[ have proposed a 
platform, based on the linearized radiation-pressure in¬ 
teraction, to implement general Gaussian operations [25l — 
|27 | between multiple mechanical oscillators. The imple¬ 
mentation of such a platform would represent a first step 
towards the realization of the circuit model of universal 
quantum computation over continuous variables [28j, . 

However, a valid alternative approach to quantum com¬ 
putation is constituted by the so called measurement- 
based model 0. The latter allows to perform general 
processing of quantum information over continuous vari¬ 
ables [31] provided a massively entangled state — dubbed 
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cluster state — is used as a resource and additional mea¬ 
surements are locally performed over its constituents. 
Despite the limitations of finite squeezing (3^, [33j, this 
approach has been proven to be fault tolerant [341 ] and, 
as a matter of fact, much effort has been devoted towards 
the generation of cluster states of light, culminating with 
the experimental realization of states composed ofa high 
number of modes (up to 10,000) both in the time [35] and 
frequency [3^, ^7i] domain. However, specific schemes for 
the generation of cluster states involving massive degrees 
of freedoms, rather than radiative ones, are still lacking 
despite some theoretical framework generically suit¬ 
able for their implementation has been proposed [381 . 1391 ]. 
The main advantage of this type of cluster states is that, 
being hosted in stationary or solid-state based architec¬ 
tures, they offer a promising path towards integrated and 
scalable quantum technologies. 

In order to bridge this gap, the aim of the present work 
is to introduce a scheme to generate continuous-variable 
cluster states of mechanical oscillators. We propose a 
scheme for generating graph states [40] of arbitrary size 
and shape (including cluster states) whose nodes are em¬ 
bodied by the mechanical modes of an optomechanical 
system. The graph states are obtained by properly en¬ 
gineering both the Hamiltonian and the dissipative dy¬ 
namics of the radiation degrees of freedom. Specifically, 
the method we use to engineer the desired Hamiltonian is 
based on multi-tone external dri ving , adapting and gen¬ 
eralizing previous approaches j41H4fij so that the required 
sidebands could be independently excited. In order to 
drive dissipatively the system to the graph states, we 
use a theoretical framework — introduced in Ref. |47| - 
that adapts quantum dissipation engineering to Gaussian 
continuous-variable systems. The merit of our scheme is 
that one can generate arbitrary graph states only by driv¬ 
ing the optomechanical system with a sequence of tunable 
pulses. 

The paper is organized as follows. In Sec. |H]we intro- 
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duce the system under consideration and derive the tun¬ 
able linearized Hamiltonian that will be instrumental for 
the graph-state generation. In Sec. Mil we will introduce 
the state generation protocol for the case of a generic 
graph state and in the absence of mechanical noise. We 
illustrate the action of the protocol via considering two 
specific examples. The detrimental effect of mechanical 
noise will be considered in detail in Sec. IIV1 confirming 
the robustness of the present protocol for low noise. A 
brief discussion about the experimental feasibility of our 
scheme is given in Sec. [V] Finally, Sec. IVII will close the 
paper with some concluding remarks. 


II. TUNABLE LINEARIZED HAMILTONIAN 
OF THE SYSTEM 

Consider an optomechanical system consisting of N 
non-interacting mechanical resonators and a single-mode 
optical cavity driven by M classical laser fields [M -tone 
drive, see Fig. 0(a)]. This configuration extends and 
adapts to our purposes some approaches already consid¬ 
ered in the literature in order to generate single-mode 
squeezing 1411 and entanglement among two mechanical 
oscillators [42 - l46l |. In this section, we derive the fam¬ 
ily of tunable Hamiltonians that we will use later on to 
generate our target states. 


to the j th mechanical resonator mode. The driving laser 
fields e(t) are given by 

M 

e(t) = Y, efce-^V** (2) 

fc=l 

where u>k, e*, and &k are respectively the frequency, in¬ 
tensity and phase of the k th laser field. 

The equation of motion for the optical and mechanical 
modes are written as: 

N 

a = — iu> c a — ia 9j(^] + bj) — ie(t) 

i=i 

--a-v^Oin (3) 

bj = —ittjbj — igja^a (4) 

where ai n is the input noise operator for the cavity mode 
with decay rate n. In Eq. 0 we have omitted the terms 
describing the mechanical noise, assuming that the cav¬ 
ity damping process dominates all other noisy processes. 
This requirement will be relaxed in Sec. ED where we will 
consider the robustness of our scheme against mechanical 
noise. 

Now, we linearise the equations of motion by consid¬ 
ering the following: 



FIG. 1. (Color online) (a) : An optomechanical system 
consisting of one optical cavity mode a coupled to N non¬ 
interacting mechanical resonators 6i,..., 6jv. The cavity dis¬ 
sipates with a damping rate k, and it is driven by M classical 
laser fields, (b) : The state of the mechanical resonators can 
be prepared in different graph state geometries, e.g., from left 
to right, a linear, a dual-rail, and a generic graph state. 


The Hamiltonian of the system described in Fig. [T] is 
given by fl], 0 


N 

H = LOcO 1 a+Y^ 
j =i 


+ gjata(bl+ bj) +e(t)a^+e* (t)c 


( 1 ) 

where a and u> c (bj and flj) are the annihilation operator 
and frequency of the cavity (j th mechanical-resonator) 
mode respectively, whereas gj is the one photon coupling 


gj<€. f - j 

(5) 

a—» a + a 

(6) 

bj^f bj + Pj 

(7) 

where, in Eqs. (0 and 0, we split the optical (mechan¬ 
ical) field into classical part a (fij) and a quantum noise 

term a (bj). Neglecting the second-order terms 
bj , the equations of motion become: 

in a and 

N 


a « -~iu] c a 9j (Pj + Pj) 

j = 1 


N 


~ia'Y J 9j{b) + bj) - -<a - \fKa in 

7 = 1 

(8) 

J A - 

N 

K 

-a (9) 

a « — iuj c a — ia 9j(Pj + Pj) — *e(t) — 

j=l 

bj ps —iVLjbj — igj(aa 1 + a*a) 

(10) 

Pj ~ ~ *£b'H 2 

(11) 


and the linearised Hamiltonian reads: 

N N 

H = w' c a 1l a+'Y, fybjbj + 'Y / 9j{ a a^ + at*a)(tfj+bj) ( 12 ) 

j=i 3=1 

with w' = w c + J2jLi9j(Pj + Pj)- For \Pj\gj < w c we 
have oj' c « ui c . 
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From Eq. ©, the parameter a is approximately given 

by 


« = E 


-iek 
k /2 + i(ui c 


Wfc) 




(13) 


By defining a k e i<t>k := K / 2+i l^ IJJk ) ^ k with a k > 0 and 
4>k £ IR, the linearised Hamiltonian becomes: 


N 


H = uj c a)a + fljb^bj 
i=i 

N M 

+ EE ®k9j& fc + fy) “I - H.c. 

j=i fe=i 


.(14) 


In order to eliminate the explicit time dependence of 
the Hamiltonian we first go to the interaction picture 

N M 

H = a) EE akgje^e^ c k ^ ^e j + e j bj'j 

3=1 fe= 1 

+H.c. (15) 


In general, the dissipation into an inaccessible environ¬ 
ment causes the loss of quantum coherence, driving the 
system of interest into a state void of genuine quantum 
features. However, dissipation can indeed turn into a re¬ 
source that does stabilise quantum coherence, provided 
the dissipative dynamics is properly engineered [48i - f52j . 
In the context of continuous-variable systems, dissipation 
engineering has been considered for the purpose of entan¬ 
glem ent generation, both at a general theoretical level 
55|] and in experiments involving atomic ensem- 
-I58j ■ Concerning the specific setting of optorne- 
chanics, recent theoretical proposals have shown how dis¬ 
sipation engineering can allow to achieve strong steady- 
state squeezing [4l| and two-mode entanglement in three- 
and four-mode systems (42l-l4?it Ififll. |60|. Here, we gener¬ 
alize these approaches and propose a scheme to achieve 
an arbitrary graph state at the steady state. 

The dynamics of the system is described by the master 
equation 


dp 

At 


-i[H 1 p\ + C 


(18) 


We then consider for each mechanical resonator j, two 
laser drives with frequencies uj± = w c ± f lj , amplitudes 
a±, and phases (fr). Assuming that the mechanical fre¬ 
quencies do not overlap and that we are in the weak cou¬ 
pling regime 


with p is the system’s density matrix, H is given by 
Eq. (fI71) . and C accounts for the dissipation processes. 
Assuming that the cavity damping process dominates all 
other dissipative processes, C reads as follows: 


afgj « , (16) 

we can invoke the rotating wave approximation, finally 
obtaining the following Hamiltonian: 

N 

H = a f E 9j (a+e^ + b] + ajeb + H.c. (17) 

3 = 1 

The main feature of the Hamiltonian above is that it 
includes both “beam-splitter” and “two-mode squeezing” 
interactions that can be independently tuned for each me¬ 
chanical oscillator. In particular, we will see in the next 
Section that, by controlling the amplitude and relative 
phase of each driving laser, one can conveniently tune 
the parameters a± and to generate the desired graph 
state. 


III. GRAPH-STATE GENERATION PROTOCOL 

As said, our aim is to devise a scheme to generate arbi¬ 
trary graph states of the mechanical degrees of freedom 
for the system described in Sec. [TT] [see Fig. H}(b)]. We 
will show now how this goal can be achieved by exploit¬ 
ing as a resource the dissipation of the radiation mode, 
together with a suitably tuned sequence of laser pulses. 
For the sake of clarity, in this Section we will neglect the 
dissipation of the mechanical degrees of freedom, which 
will be taken into account in Sec. hyi 


C = n(apa^ — —a^ap — —pa) a). (19) 


One can show that the system above admits a unique 
pure Gaussian steady-state only if just a single mechan¬ 
ical mode is involved (see Appendix [A| • This obstacle 
can be circumvented, in the case of many mechanical res¬ 
onators, using the mechanism of Hamiltonian switching 
proposed by Li et al. [6l| and generalised by Ikeda and 
Yamamoto [47|- The key point of this approach is that 
only one suitable collective mode is in fact coupled at 
each switching step. This is achieved by dividing the dy¬ 
namics in as many steps as the number N of (mechanical) 
modes involved. At each step, the interaction Hamilto¬ 
nian has to be properly engineered in order for the op¬ 
tical dissipation to cool a specific collective mechanical 
mode. Once the steady state is achieved, another collec¬ 
tive mode is cooled by switching to a new Hamiltonian 
(with the remaining collective modes uncoupled). The 
main difficulty of this scheme is to design a system that 
sustains the set of Hamiltonians required to generate a 
certain desired state, in our case an arbitrary Gaussian 
graph state. 

More specifically, consider a generic quadratic uni¬ 
tary transformation U l25l-f27f acting as c = Ub 
on the mechanical modes b = (bi ,..., b]\r) T , where 
c = (ci,..., cn) t defines a set of collective modes. As 
said, we perform an 7V-step transformation, where in the 
fc th step we set the driving laser’s amplitudes and phases 
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Step 

«i 

Oi 2 

«3 

« 4 


^2 

03 

04 

1 

\/ 2(5+\/5) 

5 

\/5+2v / 5 

5 

\/ 5—2\/5 

5 

\/5+2V5 

5 

3ti-/2 

TV 

7t/2 

0 

2 

\J 5+2\/5 

5 

y/s+2V5 

5 

\J 2(5—\/5) 

5 

-\/5-2V5 

5 

7V 

3tt/2 

TV 

tt/2 

3 

\/5-2\/5 

\/2(5 —\/5) 

\/ 5+2\/5 

V / 5+2\/5 

71-/2 

TV 

3n/2 

TV 

4 

5 

\/5+2V5 

5 

V2(5-V5) 

5 

\/. 5-2n/6 

5 

0 

7t/2 

TV 

3tt/2 


TABLE I. Assignement of the parameters in Eq. (1171) for each of the four steps of the Hamiltonian switching procedure that 
allows to generate a four-mode linear graph state. 
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«3 
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0i 
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3 

4 

5+x/5 

10 

1 

-5+V5 

10 

1 

1 

V5 

5+VS 

10 

1 

y/5 

-5+V5 

10 

-5+x/5 

10 

1 

V5 

5+x/5 

10 

1 

1 

Vs 

- 5+VS 

10 

1 

%/5 

5+V5 

10 

3tt/2 

TV 

tv/2 

TV 

TV 

3tt/2 

71 

71/2 

tv/2 

TV 

3tt/2 

TV 

TV 

tv/2 

TV 

3n/2 


TABLE II. Assignement of the parameters in Eq. (El for each of the four steps of the Hamiltonian switching procedure that 
allows to generate a four-mode square graph. 


in Eq. m as follows: 

a J = ~~\Ukj\ (20) 

9j 

aj~ = raj (21) 

4>j = = arg {U kj ) (22) 

where /3 > 0 and 0 < r < 1. With this settings Eq. 1171) 
becomes: 

H = H^ =pa\c k +rcl) + H.c. (23) 

One can show (see Appendix [Aj that the above Hamilto¬ 

nian, with the help of the cavity dissipation in Eq. m, 
generates at the steady state a single-mode squeezed 
state for the collective mode Ck- After applying all 
the N steps of the transformation, we obtain N single¬ 
mode squeezed state relative to the modes ci,...,cjv- 
In terms of the mechanical modes b\, ..., bisr, the latter 
corresponds to a Gaussian state with zero mean. Let 
us introduce the canonical position and momentum op¬ 
erators for each mechanical mode qj = (fet + bj)/y/ 2 
and pj = i(bj — bj)/y/ 2, respectively, the vector R = 
(qi ,..., qN,Pi, ■ ■ ■ ,Pn), and the covariance matrix V 
whose elements are [V]ki = (RkRi + RiRk)/ 2. Then 
the covariance matrix of the mechanical Gaussian state 
reads: 



e -2£-^NxN q NxN 
QiVxiV e 2 t-^NxN 


S 


(24) 


(l24l) has a level of squeezing of 10 log 10 [e 2 ^]dB. Ikeda 
and Yamamoto [47} showed that the unitary transfor¬ 
mation U that generates an arbitrary Gaussian graph 
state is obtained from the polar decomposition of the 
matrix — (*1 jv + A) = R.U where A is the adjacency ma¬ 
trix describing the graph state that one wants to generate 
m , and R is N x N real matrix. 


A. Examples 

In the remaining part of this section, we give ex¬ 
plicit examples demonstrating how to generate a linear 
and a square four-mode graph states using the foregoing 
scheme. We have chosen four modes only for the sake 
of clarity but, in principle, the size and shape of graph 
state that can be engineered is arbitrary. For simplicity 
we set gj = g for all j = 1,..., 4, and all the amplitudes 
are given in units of (3/g. 

Let us first consider a four-mode equally weighted lin¬ 
ear graph state. This state is a basic building block in 
measurement-based quantum computation since it allows 
to implement an arbitrary single-mode Bogoliubov trans¬ 
formation via local measurements [g4|,(65j- The adjacency 
matrix is given by: 

/ 0 1 0 0 \ 

A= 1010 

0 10 1 

\0 0 1 0 / 


where £ = tanh -1 r and S is the sympletic transforma¬ 
tion corresponding to the quadratic unitary transforma¬ 
tion U [25l-f27j|. The quantity e _2? is the squeezing pa¬ 
rameter and the state given by the covariance matrix 


from which, as explained above, we obtain the corre¬ 
sponding unitary matrix U via the polar decomposition 
of the matrix — (il .4 + A). The graph state is then gen¬ 
erated in four steps. In each step the parameters of the 
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driving laser fields are tuned to the values obtained from 
Eqs. ([20|) and (121 and indicated in Tab. U [aj is always 
given by Eq. Ell) and </>+ = 

As a second example, we consider a four-mode equally 
weighted square graph state. The relevance of the latter 
stems from the fact that it allows to perform two-nrode 
quantum operations and, in addition, it can be used in a 
redundant encoding scheme for error filtration [(id . The 
corresponding adjacency matrix is : 

/ 0 1 0 1 \ 

A= 1010 . 

0 10 1 

V i o i o / 

As before, the graph state is generated in four steps, 
where in each step, the parameters of the laser fields are 
set according to the values indicated in Tab. HU 

In principle, the results derived here are strictly valid 
only in the ideal case in which the steady state is reached 
at each step of the Hamiltonian switching scheme. In 
practice, the amount of time that can be devoted to each 
step is finite and one should assess the errors determined 
by the finite-time dynamics. The latter can be estimated 
using the Uhlmann fidelity between the ideal target state 
and the actual one, as a function of the switching time t s 
(ie., the time elapsing between one step and the other). 
Given two generic states p\ and p 2 , the Uhlmann fidelity 
is defined as F(pi,p 2 ) = (Tr y 7 y/piP 2 y/pi) > which sat¬ 
isfies 0 < F(pi,p 2 ) < 1 and F{p-\ , Pi) = 1 if and only 
if the states p\ and p 2 are identical [63]. The latter can 
be conveniently expressed in terms of covariance matri¬ 
ces for Gaussian states, using the expressions provided 
in Refs. [H, In Fig. [2] we plot the time evolution of 
the fidelity for a finite switching time t s and for the case 
of a four-mode linear cluster state. We can see that the 
fidelity approaches its final value F gen step after step. In 
Fig. El we report instead the final fidelity F gen as a func¬ 
tion of t s , for the same target state. We can see that, as 
expected, the fidelity of the state generated via Hamil¬ 
tonian switching and dissipation engineering approaches 
the unit value for large enough switching times t s . In 
Appendix [Bj we give the time scale to reach the steady 
state at each step of the switching scheme. In the nu¬ 
merical simulations shown in Fig. [2] and Fig. [3[ we have 
chosen /? = k/( 4V1 — r 2 ), so that the system reaches 
the steady state at each step in minimal time (see Ap¬ 
pendix [Bj and with minimum driving power [the driving 
power determines the value of otj. which are related to 
/? via Eq. EH1) ]. With this choice of /3, using Eq. EU1) 
and considering the condition gsd, the cavity decay rate 
must satisfy the following condition: 

k -C 4\/l — r 2 min / . (25) 

The latter implies that (i) the scheme here introduced 
works in the resolved-sideband regime and (it) the higher 


the squeezing the deeper into the good cavity regime one 
has to operate. 



FIG. 2. (Color online) Time evolution of the fidelity for 
a four-modes linear graph state with fixed switching time 
t s = 20kg 1 . We have set for the j th oscillator the frequency 
Qj/ 2tt = j MHz. 



FIG. 3. (Color online) Fidelity of the steady state of the 
four-mode linear graph as a function of evolution time per 
switching step. The fidelity shown in this plot is that of the 
steady state of the mechanical oscillators after applying all 
the switching steps. We have set for the j th oscillator the 
frequency Qj/2n = j MHz. 


IV. ROBUSTNESS OF THE GRAPH-STATE 
GENERATION AGAINST MECHANICAL NOISE 

We have hitherto showed how to generate the tar¬ 
get state without considering any mechanical dissipation. 
We now relax this condition and consider the mechanical 
resonators in contact with a thermal bath, addressing the 
detrimental effect of the latter to the scheme introduced 
above. 

In the master equation (usd, we include the following 
additional decoherence channels: 

N ( , 1 , i \ 

A = ^2lj( n j + !) ( b iP b ] - 2 b ) b jP~ 2 pb fo ) ( 26 ) 

.7=1 V ' 

N ( 1 1 \ 

A 7i n 3 ( b )p b o - 2 bjb j p ~ 2 pbjb 

. 7=1 ^ ' 
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Temperature (Kelvin) 


FIG. 4. (Color online) Contour plots of the fidelity between the target state and the state generated using the protocol described 
in the text. The fidelity is shown as a function of the temperature of the mechanical bath (horizontal axis) and the mechanical 
damping rate (vertical axis) for different target graph states and levels of squeezing. Each mechanical oscillator has a frequency 
flj/2n = j MHz (j = 1,..., N, with IV = 2, 4, 8). The solid white lines correspond to a fidelity of 0.99, the dashed lines to 0.90, 
and the dotted lines to 0.80. Each data point is taken for an optimal evolution time (see text). 


with 7 j <C fij , where jj and rij = ^exp — 1J 

are the damping rate and the mean phonon number of 
the j th mechanical bath respectively, with Tj being the 
bath temperatures. Notice that we are assuming high- 
Q mechanical oscillators, hence we have considered the 
Markovian master equation above |4l j-l44l] rather then the 
full non-Markovian Brownian noise [70J. 

The thermal baths interact with the system during the 
preparation of the graph states, which results in a devi¬ 
ation of the final state from the target one. To quan¬ 
tify this deviation, we use again the Uhlmann fidelity. 
To see the effect of the thermal noise on our results of 
Sec. mi we calculate the fidelity as a function of the 
mechanical damping rates and the bath temperatures. 
In our numerical simulations, we took the latter equal 
for each mechanical oscillator (7 = 71 = • • • = 7 n , 
T = T\ = ... = Tjsi) and set mechanical modes fre¬ 
quencies as Q.j/2-K = j MHz (j = 1 Figure Q] 

shows the fidelity for 2, 4 and 8-modes graph states as 
a function of "//k and T. Each point of our simula¬ 
tion was obtained by searching for the optimal evolu¬ 
tion time of the system in order to have maximum fi¬ 
delity: very short evolution times do not allow to reach 
the target graph state, while very long times lead to dom¬ 
inant decoherence effects caused by the thermal bath. 
In other words, in the presence of mechanical noise, 
the best fidelity is not achieved at the steady state but 
at an intermediate evolution time typically larger than 


[ft (f - ~ I3' 2 (l-r 2 )) ]- 1 (see Appendix 0) so 

that the desired target state can be achieved — and small 
compared to N times the inverse of the thermal decoher¬ 
ence rate i.e., smaller then mini< 7 '<jv 1 so that the 

mechanical noise is not too detrimental during the N 
switching steps. 

One might wonder whether the fact that we are no 
longer considering the steady state constitutes an issue. 
This is in fact not the case in the context of cluster state 
quantum computation, where it is not necessary to build 
the entire cluster ahead of the computational task to per¬ 
form, but rather it is possible to build and consume the 
cluster during the computation itself (see, e.g ., Refs. [zi^ 

H). 

In Fig.[4]we focused on three paradigmatic graph states 
and squeezing levels. The first row shows results for the 
simplest graph, consisting of two modes only. Up to lo¬ 
cal symplectic operations that do not change its entan¬ 
glement, the latter is equivalent to the usual two-mode 
squeezed state whose generation in optomechanical sys¬ 
tems have been already addressed in previous literature 
[flM Ififll . (6(| . The second row refers to a four-mode 
linear graph, which as said allows to implement a gen¬ 
eral single-mode symplectic transformation. The third 
row concerns an eight-mode dual rail graph, which al¬ 
lows instead to implement a generic two-mode symplectic 
transformation. The last two graphs thus encompass the 
necessary building blocks for universal multi-mode Gaus- 
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sian unitaries. Each column in Fig. [prefers to a specific 
squeezing level for the target state: from left to right 
we set 5dB, 12.7dB, and 21dB. The first level of squeez¬ 
ing coincides with the one of the largest optical cluster 
state reported to date [35}, whereas 12.7dB is the high¬ 
est optical squeezing experimentally achieved with one 
optical mode [74}. The third squeezing level is a theo¬ 
retical upper bound to the squeezing required to perform 
universal fault-tolerant quantum computation over con¬ 
tinuous variables (see Ref. [34}). As before, the plots of 
Fig. U were obtained by setting /3 = k/( 4-\ZT^r^) in or¬ 
der to reach the maximum fidelity in minimum time (see 
Appendix iBl for more details). As said, this implies that 
the cavity decay rate must satisfy condition (ESI). 

From Fig. U we can see that, for a given squeezing and 
number of modes of the target state, the fidelity between 
the latter and the state obtained with our protocol de¬ 
creases as the temperature T and the coupling with the 
bath 7 increase. This is to be expected since our protocol 
will end up in a steady state that is far from a pure state 
for large mechanical noise. In general, one can also see 
in Fig. [I] that the region of high fidelity shrinks from the 
left panels to the right ones. This means that, for a given 
graph, the higher the squeezing of the target the less the 
mechanical noise that can be tolerated. The same be¬ 
haviour can be see from top to bottom, implying that for 
a given target squeezing the larger the graph state the 
less the mechanical noise tolerated. 

Let us notice here that the choice of the frequencies 
used in Figs. [I] and 0 is not unique. Indeed, our proto¬ 
col works as well for any other choices of the oscillators 
frequencies, as long as they do not overlap and the rotat¬ 
ing wave approximation can therefore be applied to the 
Hamiltonian in Eq. m- 


V. EXPERIMENTAL FEASIBILITY 

In view of the results of the previous Sections, let us 
now discuss some aspects regarding the experimental fea¬ 
sibility of the state generation scheme. 

First, let us consider the attainability of the system 
Hamiltonian given in Eq. m- As mentioned, vari¬ 
ous experiments have recently succeeded in realizing the 
weak optomechanical interaction that we have considered 
here [1H3}. The main requirement that differentiates our 
scheme from the latter is that we consider, rather than 
only one mechanical oscillator, multiple oscillators with 
non-overlapping frequencies. The first implementations 
of such systems have been reported recently [751477} , thus 
providing a promising route towards the realization of 
small optomechanical networks. In addition, in quan¬ 
tum electromechanical systems, mechanically complaint 
membranes of slightly different geometry and size allow 
for the realization of mechanical oscillators with non over¬ 
lapping frequencies — as required in order for the rotat¬ 
ing wave approximation adopted in Eq. o to be valid. 
For example, the experiments described in Refs. 0. 78 


111 report a mechanical frequency of 3.6 MHz, 4 MHz, 
10.56 MHz, and 13.03 MHz respectively. This wide range 
of frequencies, obtained in compatible experimental set¬ 
ups, conveniently suits our proposal. More in general, the 
variety of mechanical frequencies realized in optomechan¬ 
ical systems (even within similar settings) suggests that 
the realization of small optomechanical networks with 
non-overlapping frequencies should be within reach. 

We should mention here that, despite the fact that our 
results are obtained for non-interacting non-degenerate 
mechanical resonators, the scheme here introduced can 
be easily extended to the case of interacting mechanical 
modes. In this case, the cavity mode could interact with 
one mechanical oscillator only (see Fig. E3). After diag- 
onalising the mechanical Hamiltonian, one obtains a set 
of non-interacting mechanical normal modes with non¬ 
overlapping frequencies, thus recovering the case con¬ 
sidered here. The non-degeneracy of the normal modes 
could in principle be enforced by controlling the coupling 
between the mechanical modes. 



FIG. 5. (Color online) An optomechanical system consist¬ 
ing of one optical cavity mode and N mechanical resonators 
modes. The cavity mode couples only to the first mechanical 
mode, and all the mechanical modes couple to each other via 
nearest-neighbour interaction. 

Given the findings of Sec. El our generation protocol 
performs better in the resolved-sideband regime and for 
high quality-factor oscillators at low temperature. High 
fidelity can then be achieved when 7 C k C H and 
T <C 1, a regime that has been extensively considered for 
sideband cooling in quantum electromechanical systems. 
For example, using experimental parameters of the order 
of the ones of Ref. [5|, the values of fidelity shown in 
Fig. E] could be achieved. 

Finally, let us briefly mention a possible readout strat¬ 
egy suitable to the setting here considered. In general, 
various approaches have been put forward in order to re¬ 
construct the state of a single mechanical oscillator in an 
optomechanical setting (see Ref. [81] for a recent review). 
The generalization of those approaches to many mechan¬ 
ical oscillators could follow the path already pursued in 
other experimental contexts, such as ions or circuit-QED 
systems. In particular, in Ref. [82] the state of a net¬ 
work of interacting oscillators is reconstructed via cou¬ 
pling only one of them to a two-level system, which plays 
the role of a probe that is eventually measured. The 
optomechanical setting here considered does not include 
any two-level system, nonetheless the scheme of Ref. [§2| 
could be adapted to the present case as well where now 
the radiation, rather than a two-level system, could act 









b) 



FIG. 6. (Color online) Fidelity as a function of the target 
squeezing between a linear graph state of N nodes and the 
state generated using the protocol described in the text. Each 
data point is taken for an optimal evolution time (see text). 
Each mechanical oscillator has a frequency D.j/2n = 11 j MHz 
(j = 1 ,,N, with N = 1,..., 6) and mechanical damping 
7/2-7T = 32 Hz. The cavity mode decays with rate k/2ty = 
0.2 MHz. The bath temperature for panel (a) is T = 15 mK 
and for panel (b) is T — 1 mK 

as a probe. In fact, a tomographic scheme along these 
lines is currently under study |83|. In particular, it is 
possible to show that such a scheme works both for the 
case of interacting and non-interacting mechanical oscil¬ 
lators in an optomechanical setting (ie., both in the case 
of Fig [T]and Fig [5]). 

VI. CONCLUSIONS 

We have shown how to prepare arbitrary graph states 
of mechanical oscillators in the optomechanical system of 
FigJTI The preparation is achieved by properly driving 
with external fields each of the two sidebands correspond¬ 
ing to each mechanical oscillator. The target graph state 
is generated by switching on and off a sequence of lin¬ 
earized Hamiltonians, which in turn can be achieved by 
changing the intensities and phases of the driving fields. 
At each switching step, the external pulses are applied 
for a sufficient time until a steady state is obtained by 
virtue of the optical losses. 

In addition, we have considered the effect of non¬ 
zero mechanical noise, confirming the robustness of the 
present scheme. For a low number of mechanical oscil¬ 
lators and moderate squeezing, our protocol appears to 
be within reach of current technology since it requires 


linearized radiation-pressure interaction, resolved side¬ 
band, and low mechanical noise — a regime that has 
been already achieved for the case of a single mechanical 
oscillato r [5 | and approached for the case of few oscilla¬ 
tors [751477]. In addition, experiments have been recently 
reported Mm where single-mode mechanical squeez¬ 
ing is achieved using the approach proposed in Ref. j4l| . 
This is, in turn, very promising for the realization of the 
state generation protocol here presented, given that the 
latter extends and adapts the approach of Ref. jrjj to 
generic multi-mode graph states. This work therefore 
identifies a promising path towards cluster states genera¬ 
tion in mechanical systems, thus representing a first step 
towards measurement-based computation over continu¬ 
ous variables in a solid-state platform, rather then in the 
common optical setting. 
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Appendix A: Existence and uniqueness of Gaussian 
steady states 

Consider an open system of N modes with Hamiltonian 
H and dissipation channels L 2l ■ ■ ■, iif • The system 
evolves in time according to the master equation: 

M 

p = —i[H, p]+J2 (- L k pLl - 1/2 L\L kP - 1/2 pL\L k ) . 

k 

(Al) 

We define H = 1/2 R T GR and L = 

(Li, L 2 ,..., Lm) t = CR with G=G t G IR 2Arx2iV , c g 
<L Mx2N and R = (qi, q 2 ,. ■., qw, Pi, P2, ■ ■ ■, Pn) T ■ 

For a Gaussian state of covariance matrix V and mean 
value (R ), the master equation (1A1I) is equivalent to |54| : 

5^ = AV + VA T + B (A2) 

at 

^§- = A(R) (A3) 

where the matrices A and B are given by: 

A = E[G + 3(C t C)] (A4) 

B = E]R(G t G)E T (A5) 

and 

( qNxN qNxN 
_^NxN qNxN 
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where U iYxAr and 0 NxN are the identity and zero matri¬ 
ces respectively, and the superscript denotes the dimen¬ 
sion of the matrix. Equation (ED has a unique steady 
Gaussian state if and only if A is Hurwitz matrix i.e., 
real part of each eigenvalue is negative. 

For the system described in Sec. HT1 given by the vector 
R = (qi, ...,q N , <?cavity,Pi, ■ ■ ■ ,PAr,P cav ity) T , the Hamil¬ 
tonian is given by equation m, and the matrices G and 
C are given in the block matrix form: 


G = 


/ qNxN 

(A Nxl ) T 

q NxN 

V (c Nxl f 


c = 



^Afxl q NxN 

0 V lxN 

xNjT qN x N 

0 (B Nx1 ) t 

1 , O 1 ^ , i) 


qN x 1 
0 

gNxl 

0 


\ 

A 6 ) 

> 


(A7) 


The eigenvalues of the matrix A are: 

rW © 2 -- 4 ’' B+p - c ■ 

“ T' " 2(jv-i) times degenerate 

2 times degenerate 

It is clear that when TV > 1, the matrix A is not Hurwitz 
and therefore there is no unique Gaussian steady state 
for the system. 

For TV = 1, we have: 

x± = -i ± \l(iY + ° 2 [( a+ ) 2 -(«-) 2 ] 

where a ± are the amplitudes of the two driving lasers. 

There exists a unique steady state if and only if a + < 
a ~, and it is obtained after a time of the order of r = 

(see Appendix ITU;. 


where, as before, the superscript denotes the dimension 
of the corresponding matrix, and the matrix elements are 
found to be as follows: 

Af x 1 = gj (at cos 0 + + aj cos j>j) (A 8 ) 

Bf xl = gj(-a1j cos 0+ + a~ cos 4>~) (A9) 

Cf xl = gj (oA sin 0+ + aj sin </>“) (A10) 

T>j XN = ^(a+sin^t - a j sin^“) (All) 

I 


Appendix B: Time scale to approach the steady 
state 

The solution of the differential equation (IA2I) is: 

V(t) = e A ^- to ) F 0 e AT(t_to) + [ ds e j 4 (t - s) He AT(t " s) 

Jto 

(Bl) 

where Vo is the state of the system at time t 0 . Diago¬ 
nalizing the matrix A = PDF -1 , with P given by the 
eigenvectors of A and D = Diag(Ai, A 2 ,...) by its eigen¬ 
values Aj, Eq. (ED becomes: 


V{t) = Pe D{t - to) P~ 1 V 0 P T 1 e D (t~ t o)p T + f d s Pe D{ - t - s) P~ 1 BP T 1 e D ( t -^)p T (B2) 

Jt 0 


At the k th step of the switching scheme described in 
Sec. EH we have the Hamiltonian (BSD and dissipa- 
tor ©. Hence the matrix A has eigenvalues: 


A ± = -J±y(0 2 -/32(l- r 2 ) (B3) 

with negative real part, implying a unique steady state. 
The speed at which the system approaches its steady 
state depends on how much negative is the real part of 
the eigenvalues A 111 : the more negative the faster the ap¬ 
proach. Defining r as the time scale to reach the steady 
state, we can estimate it as: 


1 

(ftA+ 


(B4) 


and the condition of a maximum speed for approaching 
the steady state is given by: 

^ < /3\/l — r 2 . (B5) 

Under this condition, the shortest time scale to reach the 
steady state is 


_ 4 

r min - - . 


(B 6 ) 


We should mention that, in order to prepare an infinitely 
squeezed state, an infinite amount of time is required, 
and this is clear from Eq. (IB41) where the denominator 
vanishes. 
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